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Abstract 

An efficient finite-difference scheme for the 
solution of the incompressible Navier-Stokes 
equation is used to study the vortex wake of a 
rotor in hover. The solution procedure uses a 
vorticity-stream function formulation and incor- 
porates an asymptotic far-field boundary condition 
enabling the size of the computational domain to 
be reduced in comparison to other methods . The 
results from the present method are compared with 
experimental data obtained by smoke flow visuali- 
zation and hot-wire measurements for several rotor 
blade configurations. 

Introduction 

Vortex-dominated flow fields are frequently 
encountered in the flight operation of fixed- and 
rotary-wing aircraft. Complete prediction of the 
flow field is very difficult because of the dispa- 
rate length scales associated with the generation, 
interaction, and eventual decay of the vortices. 
Methods which assume the vortex interactions to be 
inviscid and represent the flow-field as a series 
of point or line vortices have proven very 
useful. Unfortunately, they suffer from the 
disadvantage of not being able to realistically 
model either the vortex core or the merging of 
vortices. Computational methods that solve the 
Navier-Stokes equations have demonstrated the 
ability to represent the continuous vorticity 
distribution that exists within the vortex core as 
well as the merging characteristics of multiple 
vortices. However, these methods suffer from the 
disadvantage of requiring large computer times and 
storage resources along with restrictions to low 
Reynolds number flow. Recently, some progress has 
been made towards an efficient solution procedure 
for solving the incompressible Navier-Stokes 
equations^' 2 ' 3 . The procedure uses a vorticity - 
vector velocity potential formulation and incorpo- 
rates an asymptotic far field boundary condition 


enabling the size of the computational domain to 
be reduced in comparison to other methods. The 

procedure has been applied successfully in the 
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numerical solution of interacting vortex rings 
and the two-dimensional roll up of vortex wakes 3 . 

The present work is an application of the 
above procedure to the analysis of vortex wake 
merging and decay for a rotor in hover. The 
development in time of the two-dimensional wakes 
calculated in reference 3 represents the stream- 
wise development of the three-dimensional wake and 
relies on the assumption that streamwise gradients 
are small. By making an analogous assumption in 
the case of a rotor in hover that the circumfer- 
ential gradients are small, the streamwise 
development of the vortex wake can be reduced to a 
time -dependent axisymmetric calculation. Results 
from the continuous vorticity field method can be 
used to improve the free wake models currently 
used in prediction methods for helicopter rotor 

A C 

aerodynamic loading ' ' . These methods usually 
employ a large number of straight-line vortex 
segments to simulate the shed vorticity in the 
wake and always contain some empiricism in deter- 
mining the details of vortex merging as occurs in 
the roll up of the trailing vorticity in the free 
wake . 

A description of the governing equations and 
numerical procedure along with some computational 
results are presented below. Results from the 
present method are then compared with the experi- 
mental data obtained by smoke flow visualization 

7 

and hot-wire measurements . 

Governing Equations 

The flow is assumed to be governed by the 
incompressible laminar Navier-Stokes equations. 
Using a fixed inertial frame with cylindrical 
coordinates r,0, and z (r and 0 define the plane 
of the rotor rotation about the vertical axis z) 
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(9) 


with corresponding inertial velocities u, v, and 
w, the flow field induced by a rotor in hover is 
unsteady and three-dimensional. The flow field 
appears steady and three-dimensional after 
transformation to a non-inertial frame rotating 
with the rotor blade as below 

r - r 

<j> = 6 - ftt ( 1 ) 

z - z 
t - t 


The rotor blade rotates with angular rate ft about 
the z axis and corresponds to <J> = constant. The 
three-dimensional flow field can be reduced to an 
unsteady axisymmetric flow field by assuming ( 1 ) 
the circumferential velocity component and 
pressure gradient are zero and (2) the rate of 
change of the velocity gradients in the circum- 
ferential direction are small as below 


v 


_! 3£ ^ jl 

r 3<J> 2 

r 



!_ jljff 

r 2 3» 2 


= 0 


( 2 ) 


This assumption is a straight-forward extension of 
a practice commonly used in fixed-wing analyses to 
reduce the steady three-dimensional vortex wake 
roll-up calculation to an unsteady two-dimensional 
calculation. The development in time of the two- 
dimensional (in the case of fixed wing) or axi- 
symmetric (in the case of a rotor in hover) flow 
field represents the streamwise development of the 
three-dimensional flow. The resulting axisym- 
metric continuity and momentum equations are given 
below 



The vorticity is defined as 

- _3u 3w 

C 3z ~ 3r 


and the velocities are obtained as 
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The assumption that v = 0 essentially limits 
the analysis to locations outside the boundary 
layer near the rotor blade. Note that the no-slip 
condition on the rotor blade corresponds 
to v - +ftr. Equations (8) and (9) correspond to 
an initial value problem and hence the vorticity 
must be initialized at <|> - 0 , corresponding 
physically to a location a small distance 
circumferentially downstream of the rotor blade 
trailing edge. The initial vorticity distribution 
must be obtained by other means such as experiment 
or numerical lifting line theory. The rotor blade 
loading is used to determine the initial vorticity 
and a comparison of rotor blade loading obtained 
from experiment 7 with that obtained from lifting 
surface theory 4 for a two-bladed rotor 
( radius/chord - 6.0) is shown in figure 1. 


Computational Method 

The actual helical pattern of the shed vorti- 
city is approximated by a series of axisymmetric 
vorticity distributions. Their respective 
vertical locations correspond to that which would 
be obtained by passing a plane <j> =constant 
through the rotor downwash field as shown in 
figure 2. 
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where the axisymmetric Laplacian operator is 
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Note the above equations are time-like in the 
variable -<f>/ft representing circumferential 
development of the flow. The above equations can 
be efficiently solved in vorticity-stream function 
form as below 

~ a T$ + h (uc) + i <wc) = v [v2? -V (8) 

Y r 


The boundary conditions used in the far -.field 
away from the downwash field directly below the 
rotor are 

= 0(d ^) as d » (r^+z^)^*- 00 (12) 

£ = 0(e ^ ) as d = (r 2 +z 2 )^2-»- oo (X 3 ) 

corresponding to the vorticity distribution being 
of finite extent and having a distribution anala- 
gous to a classical Lamb vortex. As pointed out 
in references 1 and 2, the condition (13) is much 
stronger than (12) and a substantially more 
accurate solution can be obtained by imposing (13) 
and solving (8-9) rather than imposing (12) and 
solving (3-6) on any bounded domain. 

The vorticity transport equation is solved on 
the bounded domain D as shown in figure 2 subject 
to the boundary condition Z, ~ 0 on the outer boun- 
dary. The wake structure outside of the domain D 
has reached an asympotic state and could be 
approximated as a series of Lamb vortices. The 
movement and diffusion of these vortices could 
then be calculated by an asymptotic analysis as in 
reference 2 and their influence on the velocity 
field in D calculated exactly. For the present 
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Results and Discussion 


results, however, the vorticity distribution is 
contained entirely within the bounded domain D. 
The vorticity is initialized at $ « 0 from the 
prescribed rotor loading distribution and at 
intervals in time corresponding to the circumfer- 
ential distance from one blade passage to another, 
the initial distribution of vorticity is super- 
imposed on the existing vorticity distribution. 
The interval in time between blade passages corre- 
sponds to <J> « tt for a two-bladed rotor. The 
vertical separation of the superimposed vorticity 
distribution is obtained from experiment as dis- 
cussed subsequently. Results have been obtained 
for up to 8 superpositions of initial vorticity 
corresponding to 4 rotor revolutions for a two- 
bladed rotor » 

The velocity due to the vorticity field in D 
is calculated by using a fast Poisson solver to 
solve equation (9) for ij> and differentiating as in 
equation (11). The exact integral equation 
relating velocity to vorticity (usually referred 
to as the Biot-Savart relationship) could be inte- 
grated numerically to obtain the stream function 
in the computational domain. However, numerical 
integration of this equation over the entire flow 
field requires at least 0(N 4 ) operations where N 
is a typical number of computational grid points 
along either r or z. Numerically integrating to 

find only the boundary values requires O(N^) 
operations. Using a far-field representation for 
the boundary values of as derived in references 
1 and 2 requires only 0(N 2 ) operations. The far- 
field expansion uses moments of the vorticity 
distribution as given below 

(|) 2 { <er' 2 > +| <~) <cr’ 2 z’> 

3 z 2 2 2 

+ — [5 (-) -1] <Cr' z* > 

2d 2 d 

+ [f <f> 2 -1] <tr ,4 >]} + 0 (d“ 4 ). (14) 

where the moments. of vorticity are defined as 

OO 00 

<Cr ,m z ,n > = / J ?(r',z',it>) r ,m z ,n dr' dz' (15) 

— OO 

2 

Note the second moment of vorticity <C r* > is 
known to be conserved and hence can be used to 
measure the error of the numerical solution. 
Since evaluation of using the fast Poisson 
solver requires only O(NlnN) operations, the 
resulting procedure represents a very efficient 
method of solving the incompressible Navier-Stokes 
equations . 


Results using the procedure applied to a 
representative loading distribution shown in 
figure 3 ^re discussed below. The corresponding 
initial vorticity distribution was generated by 
summing a series of individual Lamb vortices along 
the computational grid as in reference 3. The 
strengths of the Lamb vortices are determined by 
the change in load distribution between grid 
points with the core radius of each vortex pre- 
scribed to be a fixed percentage (in this case, 
2*5 percent) of the rotor radius R. 

In figure 4, the distribution of vorticity at 
several circumferential positions ( — 4>/2n > down- 

stream of the rotor blade are shown for a Reynolds 
number based on H and rotor radius R of 40,000. 
The roll up and inboard movement of the vorticity 
occurs very rapidly and the merging is almost 
complete by one revolution. Note these calcula- 
tions were made without including the influence of 
the shed vorticity beneath the plane of rotor. 

The effect of the asymptotic far-field boun- 
dary condition procedure is shown in figure 5. 
The ip = 0 boundary condition clearly shows the 
effect of the finite grid extent. Imposing the 
far field boundary condition using ip = 0 with an 
equivalent accuracy to the far field expansion 
procedure would require a much larger grid size 
with a subsequent loss of accuracy for the 
interior finite difference calculations assuming 
the practical restraint of a fixed number of grid 
points . 

The grid size used for high Reynolds number 
flow computation is required to be very fine to be 
compatible with the wake thickness at high 
Reynolds number which, in the experimental data 7 , 
was on the order of 0.003R. The computational 
time increases as Reynolds number increases and 
wake thickness decreases. In the following, all 
the computations were done at Reynolds number 
equal to 40,000 and the wake thickness was taken 
as 0.025R as a trade-off between computer time and 
storage. The grid size was constant in both the 
radial and axial directions and was equal to 
0.0125R. The number of grid points was 100 and 
200 in the radial and axial directions, respec- 
tively, and required approximately 1 hour compu- 
tational time on a VAX 11/780 computer to predict 
nine revolutions of the rotor. 

The present method has been applied to three 
different test in which the wake geometries were 
measured by a hot-wire sensor. The first case® 
considered is a single blade rotor with aspect 
ratio equal to 4.8 operating at a collective pitch 
angle equal to 6.2 and rotational speed equal to 
1350 revolutions per minute. The experimental tip 
vortex path (figures 6 and 7) shows a good cor- 
relation with Landgrebe's empirical formula^. To 
include the influence of the shed vorticity 
beneath the plane of the rotor, as discussed 
above, a new layer of vorticity of the same 
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strength and distribution as the initial vorticity 
has been generated at the plane of the rotor at 
the time of blade passage. The axial separation 
was determined from the experimental data shown in 
figure 6. The computation has been carried out 
for nine revolutions. The vorticity distribution 
computed by using a lifting surface code 4 is 
presented in figure 8 as well as in Table I of the 
Appendix. Comparisons of experiment and lifting 
surface calculations 4 for the rotor loading for 
all three rotors studied have been similar to that 
shown in figure 1 and, for all the calculations 
presented, the lifting surface calculations have 
been used to determine the rotor loading. 

The vorticity at the end of each blade 
passage was about 60% of its original value and 
thus the mutual interaction has been dominated by 
the new vorticity. The first concentrated 
vorticity has been pushed out more in the radial 
direction while the second one moves inboard 
slightly. As new distributions of vorticity were 
generated due to blade passage, the first 
concentrated vorticity moved outward more and more 
as shown in figure 9. However, the concentrated 
vorticity nearest the plane of the rotor did not 
move inboard enough to match the test data which 
were obtained at high Reynolds number condition. 
The Reynolds numbers of the experiments are .5 X 
10 6 and .75 X 10 6 for the single- and the two- 
bladed rotors, respectively. The predicted wake 
geometry did not agree with the measurements but 
the general trend was encouraging. The source of 
disagreement may be attributed to either the low 
Reynolds number or the time-dependent nature of 
the present calculations which have not yet 
reached a steady state condition. It is believed, 
however, that the low Reynolds number in the 
calculations are the predominant source of the 
disagreement. Computations approaching the high 
Reynolds number flow condition could be achieved 
by using a multi-scale technique in which a fine 
grid is used around the concentration of vorticity 
while a coarse grid is used for the general flow 
domain. This technique is currently under 
development by the present authors. 

The second and third cases considered are 
two-bladed rotors, one with aspect ratio equal to 
6.0 ( FXC blade) and the other with an aspect ratio 
equal to 13.7 (1/7-Scale UH-1 blade). The vorti- 
city distribution of each case is shown in figures 
10 and 11 as well as in Table II and III, 
respectively, of the Appendix. 

Smoke visualization photographs showing the 
movement and interaction of the tip vortex cores 
are shown in figures 12 and 13 for the AR = 13.7 
blade. The motion and decay of the shed vorticity 
for each case are shown in figures 14 and 15 and 
show similar trends to that discussed previously. 


Concluding Remarks 

An efficient finite-difference scheme for the 
solution of the incompressible Navier -Stokes 
equations has been used to study the vortex wake 
of a rotor in hover. The present method has been 
applied to three different tests in which the wake 
geometries were measured by the hot-wire sensor. 
Comparisons show that the computed concentrated 
vorticity nearest the plane of the rotor did not 
move inboard enough to match the test data, but 
the general trend of the predicted wake geometry 
was encouraging. The source of disagreement may 
be attributed to either the time -dependent nature 
of the present calculations which have not yet 
reached a steady state condition or the low 
Reynolds number. It is believed, however, that 
the low Reynolds number in the calculations are 
the predominant source of the disagreement. The 
difficulty in computing the high Reynolds number 
flow condition could be eliminated by using a 
multi-scale technique in which a fine grid is used 
around the concentration of vorticity and a coarse 
grid for general flow domain. This technique is 
currently under development. 
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Appendix 

Table X. Georgia Tech Single Blade, AR - 4.8 
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0 is collective 
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pitch angle. 
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Table III. 1/7-Scale Huey Blade, AR = 13.7 

9 c = 8°; R — 3.43 ft., ft = 1,250 rpm 
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Radial station, r/R 


Figure 1,- Comparison of experimental rotor loading 
with lifting surface theory for a two- 
bladed rotor (AR = 6,0), 



Figure 2.- Bounded domain used in rotor wake 
analysis . 
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Figure 3.- Initial load distribution used in numerical experiments 
for vorticity evolution computations . 



Figure 5.- Effect of far-field boundary condition. 


Figure 4.- Streamwise development of vorticity 
for initial load distribution shown 
in figure 3. 




Figure 8.“ Blade sectional loading for Georgia Tech 
single blade , AR = 4.8, 


Figure 6.- Comparison . of experimental tip vortex 
geometry with empirical geometry for 
AR = 4,8 single-bladed rotor. 


Georgia Tech Single Blade 
0 r/R 1 0 




Test Result (ref. 8) 
Landgrebe's Empirical 
formula (ref. 9) 


Test Results (ref. 7) 
Landgrebe's Empirical 
formula (ref. 9) 



Figure 9,- Comparison of experimental tip vortex 
path locations with present results 
for an AR = 4,8 single-bladed rotor. 



Figure 10.- Blade sectional loading for AR = 6 
two-bladed rotor. 


Figure 7.- Inboard movement of tip vortex paths for 
single and two-bladed rotors. 
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